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The  adsorption  of  simple  benzene  derivatives  composed  of  a  benzene  ring  with  NOz,  CH3,  or  NH2  func¬ 
tional  groups  on  a  semiconducting  single-wall  carbon  nanotube  is  studied  using  the  density-functional  theory 
within  the  local-density  approximation.  The  effects  of  molecular  relaxation  in  the  adsorption  process  are 
obtained,  as  well  as  the  adsorption  energies  and  equilibrium  distances  for  several  molecular  locations  and 
orientations  on  the  surface.  We  find  that  all  of  these  benzene  derivatives  are  physisorbed  mainly  through  the 
interaction  of  the  n  orbitals  of  the  benzene  ring  and  those  of  the  carbon  nanotube.  These  aromatics  do  not 
change  significantly  the  carbon  nanotube’s  electronic  structure,  and  therefore  only  small  changes  in  the  nano¬ 
tube’s  properties  are  expected.  This  suggests  that  these  benzene  derivatives  are  suitable  for  noncovalent  nano¬ 
tube  functionalization  and  molecule  immobilization  on  nanotube  surfaces. 
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I.  INTRODUCTION 

Since  their  discovery,1  carbon  nanotubes  have  stimulated 
intense  experimental  and  theoretical  interest  in  their  physics, 
chemistry,  and  materials  science.  Their  unique  structure  and 
properties  make  them  suitable  for  a  variety  of  potential 
applications.2  Single-wall  nanotubes  (SWNTs)  are  multi¬ 
functional  materials  that  have  potential  for  use  as  efficient 
gas  storage  elements  and  in  battery  devices,  as  sensors  and 
electromechanical  systems  in  nanoelectronics,  and  as  bio¬ 
compatible  agents  and  sensors  in  medicine.3 

The  interaction  of  SWNTs  with  chemical  species  is  the 
main  ingredient  for  many  of  these  applications.  For  example, 
it  has  been  found  experimentally  and  theoretically  that 
SWNT  transport  properties  can  be  changed  upon  exposure  to 
02,  H2,  N02,  NH3,  COOH,  and  to  some  benzene 
compounds.4"10  The  sensitivity  of  SWNTs  to  these  materials 
combined  with  their  inherent  characteristics  such  as  their 
small  size,  good  electrical  and  mechanical  properties,  and 
high  surface  area  have  led  many  to  envision  their  applica¬ 
tions  in  sensors  and  in  the  immobilization  of  biological  sub¬ 
stances,  and  in  their  use  as  functionalized  materials.11"16 

Aromatic  compounds  interacting  with  SWNTs  are  of  par¬ 
ticular  interest.  Noncovalent  sidewall  SWNT  functionaliza¬ 
tion  with  aromatic  organic  molecules11,12  has  attracted  in¬ 
creasing  attention.  Using  SWNTs  as  sensors  for  some 
aromatics  has  also  been  demonstrated  experimentally.9,13  In 
addition,  SWNT  functionalization  with  biological  molecules 
(such  as  proteins,  DNA,  and  nucleic  acids),  which  contain 
repetitions  of  many  aromatic  rings,  provides  potential  for 
developing  drug  delivery  systems,  molecule  immobilization, 
and  sensors.14  Therefore,  a  fundamental  theoretical  and  sys¬ 
tematic  understanding  of  SWNT  interaction  with  benzene 
derivatives,  which  appear  as  building  blocks  in  large  organic 
and  biological  systems,  is  necessary. 

Cyclic  aromatics  consist  of  one  or  more  carbon  rings.  The 
simplest  example  is  benzene  C6H6.  The  planar  arrangement 
of  its  atoms  creates  the  conditions  for  filling  the  in-plane  sp 2 
hybrid  orbitals  of  C  and  creating  a  resonant  bond  shared  by 
all  six  C  atoms.  Besides  this  aromatic  bond,  there  is  a  set  of 


7 t  orbitals  transverse  to  the  molecular  plane.  Graphene  itself, 
as  well  as  metallic  and  semiconducting  SWNTs,  has  7r  orbit¬ 
als  perpendicular  to  the  surface.17  Therefore,  the  aromatic 
molecules/SWNT  systems  can  be  viewed  as  two  interacting 
7 t  systems.  The  tt-tt  interaction  is  long  ranged,  is  relatively 
weak,  and  includes  van  der  Waals  forces.  It  leads  to  stacking, 
which  is  apparent  in  the  interlayer  bonding  in  graphite  and  in 
the  SWNT  solubilization  in  aromatic  solvents18  or  in  aro¬ 
matic  polymers.19 

Benzene  derivatives  are  obtained  by  substituting  one  or 
more  hydrogen  atoms  in  the  benzene  ring  with  a  functional 
group  such  as  hydroxil  OH,  methyl  CH3,  nitro  N02,  or 
amino  NH2.  The  simple  benzene  derivatives  considered  here 
contain  only  one  functional  group.  These  derivatives  have  a 
common  7r  structure  from  the  hexagonal  carbon  ring,  but 
they  have  distinct  properties,  such  as  permanent  dipole  mo¬ 
ments  and  different  electron  affinities  from  the  functional 
groups.  For  such  aromatics,  it  has  been  observed  experimen¬ 
tally  that  the  electronic  and  transport  properties  of  SWNTs 
can  be  modified  significantly.7,8 

In  this  work,  the  adsorption  of  three  simple  benzene  de¬ 
rivatives  on  a  semiconducting  SWNT(8,0)  is  investigated  us¬ 
ing  first-principles  calculations  based  on  the  density- 
functional  theory  (DFT).  In  increasing  order  of  electron 
affinity  of  their  functional  groups,  these  are  aniline 
C6H5NH2,  toluene  C6H5CH3,  and  nitrobenzene  Cf)H3N02. 
To  understand  the  adsorption  better,  we  compare  the  proper¬ 
ties  of  these  complexes  with  the  adsorption  of  benzene  alone, 
and  we  also  compare  them  to  the  adsorption  of  simple 
closed-shell  molecules  without  a  benzene  ring  connected  to 
the  functional  groups  NH2,  CH3,  and  N02.  These  molecules 
are  amonia  NH3  (for  NH2),  methane  CH4  (for  CH3),  and 
nitromethane  CH3N02  (for  N02),  which  are  all  stable  under 
normal  conditions.  We  also  investigate  the  adsorption  of 
these  functional  groups  by  themselves,  trace  the  difference  in 
adsorption  properties  mainly  to  whether  the  valence  shells  of 
the  adsorbates  are  partially  or  fully  filled,  and  discuss  the 
contributions  to  the  adsorption  process  from  the  tt-tt  stack¬ 
ing  and  from  the  properties  of  the  functional  groups.  We  find 
that  the  interaction  for  the  compounds  with  closed  shells  is 
always  physisorption,  and  that  in  the  cases  with  open  shells. 
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both  physisorption  and  chemisorption  are  possible. 

The  paper  is  organized  as  follows.  In  Sec.  II,  the  method 
for  the  calculations  is  described.  Section  III  gives  the  results 
from  the  calculations,  which  are  discussed  in  Sec.  IV. 

II.  METHOD 

The  results  presented  in  this  paper  are  obtained  with 
self-consistent  DFT  calculations  using  the  Vienna  Ab  Initio 
Simulation  Package  (vasp).20  This  state-of-the-art  package 
uses  a  plane-wave  basis  and  a  periodic  supercell  method.  It 
uses  either  ultrasoft  Vanderbilt  pseudopotentials  or  the 
projector-augmented  wave21  method  to  account  for  the  core 
electrons.  These  methods  require  a  relatively  small  number 
of  plane  waves  per  atom,  are  scalable  to  large  systems,  and 
are  chemically  accurate.  For  the  systems  studied  here  with 
full  valence  shells  (benzene  derivatives,  NH3,  CH4,  and 
CH3NO2),  non-spin-polarized  calculations  were  performed. 
For  the  systems  with  partially  filled  valence  shells  (NH2, 
CH3,  and  N02),  results  from  spin-polarized  calculations  are 
presented. 

DFT  is  the  most  successful  in  describing  short-ranged 
strong  forces  involving  chemical  bonds  and  large  charge 
transfer.  It  gives  imperfect  quantitative  descriptions  of  long- 
ranged,  weak  interactions  because  the  local-density  approxi¬ 
mation  (LDA)  for  the  exchange-correlation  functional  tends 
to  overestimate  the  adsorption  energy  and  the  generalized- 
gradient  approximation  (GGA)  tends  to  underestimate  it. 
DFT  does  not  include  dispersive  forces  like  the  van  der 
Waals  interaction.22,23  All-electron  quantum  chemical  calcu¬ 
lations  of  the  real  many-body  wave  function  including 
higher-order  corrections  to  the  Hartree-Fock  method  should 
give,  in  principle,  an  accurate  description  of  van  der  Waals 
binding,24  but  this  method  cannot  manage  large  molecules 
and  large  systems.  Nevertheless,  DFT  has  been  found  to  pro¬ 
vide  a  useful  model  for  large  graphitic  structures,  providing  a 
distance  between  layers  that  matches  the  experiment.  Al¬ 
though  the  formation  energy  of  such  structures  is  believed  to 
be  underestimated  by  half,  DFT  is  successful  in  the  correct 
calculation  of  the  corrugation  parallel  to  a  graphene  layer25,26 
or  along  a  nanotube.27  Another  success  is  the  qualitatively 
and  quantitatively  accurate  description  of  electronic  structure 
along  the  planes  of  graphitic  structures  and  along  the  axis  of 
carbon  nanotubes;  the  electronic  bands  in  these  directions  are 
mainly  determined  by  the  chemistry  of  sp 2  orbitals  and  do 
not  depend  significantly  on  van  der  Waals  forces.26,27  DFT 
used  together  with  semiempiric  al  potentials  for  interlayer  in¬ 
teraction  has  also  been  reported  to  improve  the  cohesion  en¬ 
ergy  of  graphite,  but  the  relative  orientation  and  distance 
between  graphitic  layers  and  the  electronic  structure  are  only 
slightly  different  from  the  DFT  results.22,23  Therefore,  due  to 
its  success  in  graphitic  compounds,  we  adopt  DFT  here  as  a 
model  for  the  study  of  the  adsorption  of  benzene  compounds 
on  SWNTs  to  gain  insight  into  their  electronic  structure. 

We  first  tested  the  DFT  method  for  the  description  of 
graphite.  The  GGA  does  not  give  a  bound  state  of  this  sys¬ 
tem.  The  LDA  gives  a  bound  system  with  an  in-plane  lattice 
parameter  a0=2.46  A  and  a  vertical  lattice  parameter  c 
=  6.71  A.  These  are  within  0.07%  and  1.3%  of  the  experi- 
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metal  values,  respectively.32  The  calculated  binding  energy  is 
28  meV/atom,  which  is  higher  than  the  energy  calculated 
previously  from  semiempirical  corrections  to  DFT.28  Com¬ 
parison  with  experimental  data  is  difficult  at  this  time  due  to 
the  varying  results  published:  measurements  by  Girifalco  and 
Lad24  give  43  meV/atom  cohesive  energy  for  graphite;  ex¬ 
periments  by  Benedict  et  al.30  show  35  meV/atom.  The  sec¬ 
ond  test  was  the  adsorption  of  benzene  on  a  graphene  sur¬ 
face.  The  GGA  yields  very  little  binding.33  The  LDA 
describes  correctly  the  adsorption  geometry  and  electronic 
structure,  which  is  the  most  important  aspect  of  this  type  of 
calculation.  The  calculated  adsorption  energy  is  almost  half 
the  experimental  thermal  desorption  energy,  which  is  consis¬ 
tent  with  other  DFT  descriptions  of  graphite.31 

Next  we  consider  the  adsorption  of  C6H6,  C6H5NH2, 
C6H5CH3,  C6H5N02,  NH3,  CH4,  and  CH3N02  on  a  semi¬ 
conducting  SWNT(8,0)  within  non-spin-polarized  LDA. 
These  systems  all  have  closed  valence  shell  structures.  Spin- 
polarized  calculations  were  performed  for  the  adsorption  of 
the  open-shell  functional  groups  NH2,  N02,  and  CH3  on  the 
SWNT(8,0). 

The  nanotube  was  fully  relaxed  before  adding  the  adsor¬ 
bates,  obtaining  an  optimized  length  of  a  unit  cell  of  4.22  A 
along  the  nanotube  axis.  Prior  to  the  calculations  for  ben¬ 
zene,  nitrobenzene,  aniline,  and  toluene,  three  unit  cells  of 
the  SWNT  along  the  tube  axis  were  also  relaxed,  obtaining 
a  total  length  of  12.67  A  along  the  nanotube  axis.  For 
the  cases  of  NH3,  CH4,  CH3N02,  NH2,  CH3,  and  N02,  a 
supercell  with  two  unit  cells  along  the  tube  axis  was 
relaxed,  obtaining  a  total  length  of  8.45  A.  The  lateral  sizes 
of  the  supercell  were  21.95  A  along  the  adsorption  direction 
x  and  17.07  A  in  the  direction  y.  In  each  case,  a  Monkhorst- 
Pack  k  grid  for  sampling  the  Brillouin  zone  was  taken  with 
seven  k  points  in  the  z  direction  and  one  k  point  in  the 
( x,y )  plane.  The  cutoff  energy  used  was  420  eV.  We  also 
checked  a  larger  Ecutojf  parameter  (£'CMfo//=460  eV)  in  some 
cases,  and  found  that  the  adsorption  energy  is  slightly  larger 
(~5  meV)  compared  to  that  for  a  smaller  Ecutofj.  The  results 
presented  here  were  obtained  with  the  smaller  Ecutoff,  which 
assured  convergence  and  was  less  demanding  on  computing 
resources. 

For  systems  studied  here,  the  energy  convergence  crite¬ 
rion  was  10-5  eV,  and  the  force  convergence  criterion  was 
0.002  eV/A.  During  the  adsorption  calculations,  the  SWNT 
was  kept  frozen  and  the  molecule  was  allowed  to  relax.  In 
order  to  test  the  effect  of  nanotube  relaxation  during  adsorp¬ 
tion,  in  some  cases  we  allowed  the  whole  system  to  relax.  In 
those  cases,  we  found  a  very  small  variation  in  energy  (a  few 
meV)  and  essentially  no  difference  in  equilibrium  orientation 
and  relative  distance  between  the  SWNT  and  the  molecule. 
For  most  molecules  the  nanotube  structure  was  essentially 
unchanged,  and  there  was  no  tendency  toward  sp3  hybridiza¬ 
tion  and  pyramidalization  of  C  atoms.  The  only  exception 
was  CH3,  for  which  relaxation  of  the  nanotube  was  also  al¬ 
lowed.  The  nanotube  surface  was  distorted,  indicating  the 
appearance  of  sp3  hybridization  in  this  case. 

For  all  molecules,  the  adsorption  energy  was  obtained 
from 
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FIG.  1.  (Color  online)  Orientations  of  C6H6  adsorbed  on 
SWNT(8,0):  (a)  hollow,  (b)  top,  and  (c)  bridge. 

Eads  —  Emol  +  Es  WNT  ~  Emoi/s\yNT ,  ( 1 ) 

where  Emol  is  the  energy  for  the  relaxed  molecule  only,  E SWNT 
is  the  energy  for  the  relaxed  pristine  SWNT,  and  EmoUSWNT  is 
the  energy  of  the  molecule/SWNT  complex. 


III.  RESULTS 

A.  Adsorption  configurations 

Several  positions  of  the  benzene  carbon  ring  relative  to 
the  SWNT  carbon  rings  were  considered  near  three  symme¬ 
try  sites  on  the  SWNT.  Figure  1  shows  these  sites — bridge, 
top,  and  hollow — and  the  adsorption  of  benzene  C6H6  on 
them.  The  adsorption  of  aniline,  nitrobenzene,  and  toluene 
was  also  studied  in  the  vicinity  of  the  same  symmetry  sites, 
where  bridge,  top,  and  hollow  refer  to  the  location  of  the 
molecular  benzene  ring  with  respect  to  the  SWNT  carbon 
rings.  Figure  2  shows  the  lowest  energy  configurations  for 
these  molecules:  near  bridge  for  aniline,  near  bridge  for  tolu¬ 
ene,  and  near  top  for  nitrobenzene.  Other  local  minima  are 
found  near  the  other  symmetry  sites,  as  listed  in  Table  I.  This 
confirms  previous  results  obtained  by  DFT-GGA 
calculations25,35  that  the  adsorption  of  aromatics  on  carbon 
nanotubes  having  a  surface  with  very  shallow  adhesion 
energy  results  in  different  local  minima  on  the  nanotube 
surface. 

We  also  performed  calculations  for  the  three  benzene  de¬ 
rivatives  in  transverse  configurations,  i.e.,  with  their  func¬ 
tional  groups  pointing  toward  the  SWNT  surface  and  with 
their  carbon  rings  situated  away  from  the  nanotube.  Different 
configurations  with  the  benzene  plane  parallel  and  perpen¬ 
dicular  to  the  nanotube  axis  were  considered.  These  corre- 
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FIG.  2.  (Color  online)  The  lowest-energy  configurations  of 
simple  benzene  derivatives  on  SWNT(8,0):  (a)  C6H5N02  near  the 
top  site,  (b)  C6H5NH2  near  the  bridge  site,  and  (c)  C6H5CH3  near 
the  bridge  site. 


spond  to  interactions  mainly  through  the  functional  groups, 
while  the  benzene  ring  has  only  an  indirect  effect  on  the 
adsorption.  For  all  of  these  cases,  we  found  that  the  adsorp¬ 
tion  energies  are  at  least  150  meV  smaller  than  when  the 
benzene  ring  is  parallel  to  the  SWNT  surface  (Figs.  1  and  2, 
and  Table  I).  This  suggests  that  the  tt-tt  stacking  is  the  most 
important  ingredient  in  the  adsorption  of  these  aromatics, 
which  was  discussed  in  other  works  [Refs.  27  and  33-35] 
among  which  the  more  critical  analysis  in  Ref.  34  showed 
that  the  tt-tt  stacking  proceeds  with  negligible  charge  trans¬ 
fer  and  hybridization.  As  we  show  in  the  next  section  for  the 
present  aromatics,  our  electronic  structure  results  are  in 
agreement  with  the  view  in  Ref  34.  The  tt-tt  stacking  found 
here  is  also  consistent  with  recent  experiments9  that  used 
electric  fields  transverse  to  the  nanotube  axis  to  attach  polar 
molecules  and  to  modulate  the  capacitance  of  the  system. 
Those  experiments  have  shown  that  even  simple  dipolar  aro¬ 
matic  molecules  (with  a  single  carbon  ring)  do  not  give  a 
significant  capacitive  response  compared  to  molecules  with¬ 
out  carbon  rings,  which  suggests  that  the  dipole  moment  of 
the  aromatics  lies  flat  on  the  surface. 

In  order  to  understand  the  adsorption  of  the  different  ben¬ 
zene  derivatives  well,  we  considered  the  adsorption  of  NH3, 
CH4,  and  CH3N02  at  the  same  symmetry  sites.  CH3N02  was 
taken  to  have  N02  toward  the  nanotube  surface.  These  are 


TABLE  I.  Benzene  derivatives  and  simple  closed-shell  molecules  on  SWNT(8,0):  adsorption  energies  and 
equilibrium  distances. 


Bridge 

Top 

Hollow 

Eads 

(meV) 

d 

(A) 

Eads 

(meV) 

d 

(A) 

Eads 

(meV) 

d 

(A) 

Benzene  (CgHg) 

260 

3.12 

249 

3.12 

151 

3.22 

Aniline  (C6H5NH2) 

302 

2.86 

283 

2.89 

241 

2.89 

Toluene  (C6H5CH3) 

264 

2.77 

259 

2.77 

172 

2.81 

Nitrobenzene  (C6H5N02) 

274 

3.03 

277 

3.00 

222 

3.05 

Ammonia  (NH3) 

75 

2.70 

84 

2.83 

113 

2.71 

Methane  (CH4) 

80 

2.83 

84 

2.72 

95 

2.79 

Nitromethane  (CH3N02),  xOy 

110 

2.93 

106 

2.98 

129 

2.88 

Nitromethane  (CH3N02),  xQz 

113 

2.98 

123 

2.91 

124 

2.96 
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FIG.  3.  (Color  online)  The  lowest-energy  configurations  for 
simple  molecules  adsorbed  on  SWNT(8,0):  (a)  CH3N02  near  the 
hollow  site,  (b)  NH3  near  the  hollow  site,  (c)  CH4  near  the  hollow 
site,  (d)  N02  near  the  top  site,  (e)  NH2  near  the  hollow  site,  and  (f) 
CH3  near  the  top  site. 

closed  valence  shell  structures  containing  the  same  func¬ 
tional  groups  as  the  simple  aromatics,  but  they  do  not  have 
the  benzene  ring.  The  adsorption  energies  and  equilibrium 
distances  are  shown  in  Table  I  and  Figs.  3(a)-3(c),  respec¬ 
tively.  In  the  lowest  energy  configurations,  NH3  and  CH4 
have  a  tripod  configuration,  where  the  N  and  C  atoms,  re¬ 
spectively,  are  on  the  nanotube  hollow  site,  and  the  H  atoms 
point  toward  the  neighboring  C  atoms  from  the  nanotube. 
For  CH3N02,  we  considered  only  the  orientations  with  N02 
toward  the  nanotube  to  obtain  a  comparison  with  C6H5N02. 
In  this  case,  we  denote  by  xOz  the  orientation  of  the  symme¬ 
try  plane  of  the  N02  in  nitromethane  parallel  to  the  nanotube 
axis  and  by  xOy  its  orientation  perpendicular  to  the  SWNT 
axis.  The  lowest  energy  configuration  is  at  the  hollow  site  in 
the  xOy  configuration  with  O  atoms  pointing  toward  two  op¬ 
posite  C  atoms  from  the  nanotube,  as  shown  in  Fig.  3(a). 

We  also  calculated  the  adsorption  energies  and  equilib¬ 
rium  distances  for  the  three  unsaturated,  open-shell  frag¬ 
ments  N02,  NH2,  and  CH3.  Since  these  radicals  have  un¬ 
paired  electrons,  spin  polarization  is  important  and  spin- 
polarized  calculations  were  done  here.  The  data  for  the  cases 
in  which  the  N  atom  points  away  from  the  nanotube  axis  are 
given  in  Table  II  for  N02  and  NH2.  Adsorption  configura¬ 
tions  are  also  possible  with  the  N  atom  pointing  toward  the 
nanotube,  but  they  have  smaller  energies  than  those  with  N 
pointing  away  from  the  tube.  Again  we  found  that  bridge, 
top,  and  hollow  are  stable  minima  for  both  xOy  and  xOz 


orientations.  For  N02  we  also  give  results  from  non-spin- 
polarized  calculations  to  illustrate  the  importance  of  spin  po¬ 
larization  in  molecules  with  an  open-shell  electronic  struc¬ 
ture.  In  Figs.  3(d)  and  3(e),  the  equilibrium  configurations 
for  NH2  and  N02  are  shown.  It  can  be  seen  that  NH2  favors 
the  structure  with  N  in  hollow,  while  N02  prefers  N  on  top 
of  the  C  with  one  O  atom  oriented  toward  the  bridge  site  and 
the  other  O  oriented  toward  the  hollow  site  of  the  nanotube. 

The  adsorption  for  CH3  was  also  calculated  with  several 
starting  configurations  on  the  bridge,  top,  and  hollow  sym¬ 
metry  positions,  with  the  C  atom  oriented  away  and  toward 
the  carbon  nanotube  with  the  inclusion  of  spin-polarization 
effects.  All  cases  converged  into  one  preferred  orientation  of 
C  sitting  on  top  of  a  C  atom  from  the  nanotube,  as  shown  in 
Fig.  3(f).  The  adsorption  energy  was  found  to  be  1.627  eV, 
and  the  equilibrium  distance  1.54  A. 

The  adsorption  energies  for  N02  and  NH2  are  relatively 
small,  indicating  physisorption.  They  are  comparable  to  the 
energies  for  the  closed-shell  compounds.  On  the  other  hand, 
the  CH3  radical  binds  much  more  strongly  to  the  nanotube 
surface  via  chemisorption. 

There  have  been  several  reports  using  DFT  for  the  adsorp¬ 
tion  on  SWNTs  of  the  small  molecules  considered  here  and 
also  for  the  adsorption  of  benzene.  For  NH3,  a  previous 
work'6  reported  binding  with  the  N  atom  toward  the  nano¬ 
tube,  while  another  work37  found  no  binding.  For  the  same 
molecule,  other  reports38,39  find  an  adsorption  energy  of 
-140  meV  on  a  (10,0)  SWNT.  For  CH4  on  a  (10,0)  SWNT, 
an  energy  of  —190  meV  and  a  distance  of  3.17  A  were 
obtained.39  For  N02,  the  adsorption  energies  obtained  using 
non-spin-polarized’740  as  well  as  spin-polarized 
calculations41,42  have  a  rather  wide  range,  between  300  and 
600  meV.  The  upper  value41  suggests  that  stronger  forces 
might  become  important.  The  chemical  adsorption  of  CH„ 
(n=l,2,3)  radicals  on  a  nanotube  has  also  been  reported 
previously,43  but  without  spin  polarization.  In  the  case  of 
benzene  on  SWNT  (8,0),  Tournus  et  al. 11  calculated  an  ad¬ 
sorption  energy  of  193  meV  and  a  distance  of  3.19  A  at  the 
bridge  site,  which  are  67  meV  lower  and  0.07  A  longer  than 
the  energy  and  distance  found  here. 

We  attribute  the  differences  between  the  results  in  these 
reports  and  ours  to  the  calculation  details.  Not  all  the 
previous  papers  include  relaxation  during  the  adsorption  pro¬ 
cess.  Here  relaxation  of  the  adsorbed  species  during  the  in¬ 
teraction  process  is  allowed.  This  explains  partially  the  dif¬ 
ference  in  adsorption  energy  for  C6H6/SWNT  between  our 
result  and  that  reported  by  Tournus  et  al ,27  In  addition,  as  a 
consequence  of  the  long-ranged  forces  governing  the  adsorp- 


TABLE  II.  Open-shell  functional  fragments  on  SWNT(8,0):  adsorption  energies  and  equilibrium  distances. 
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FIG.  4.  (Color  online)  Total 
density  of  states  for  (a) 
benzene/SWNT  (8,0),  (b)  clean/ 
SWNT(8,0),  (c)  nitrobenzene/ 
SWNT(8,0),  (d)  nitromethane/ 
SWNT  (8,0),  (e)  toluene/SWNT 
(8,0),  (f)  methane/SWNT  (8,0), 
(g)  aniline/SWNT  (8 ,0),  and  (h) 
NH3/SWNT  (8,0). 


tion,  we  found  that  the  DFT-LDA  results  for  the  current  mol¬ 
ecules  are  sensitive  to  the  value  of  the  cutoff  energy  Ecu[0fj.  A 
large  cutoff  parameter  must  be  used  to  include  many  plane 
waves  in  the  calculation,  as  described  earlier  by  Girifalco 
and  Hodat  23  for  graphite.  Using  a  small  cutoff  parameter 
may  result  in  unphysically  large  adsorption  energies.  We  per¬ 
formed  calculations  for  N02  with  Ecutofj= 420  eV  and  also 
with  Ecutofj-=  320  eV.  The  adsorption  energy  for  the  smaller 
E cutoff  was  Ealjs=5 50  meV,  which  is  in  the  range  found  in  the 
literature.41 

A  comparison  of  equilibrium  distances  for  the  four  aro¬ 
matics  studied  here  shows  that  the  distance  between  benzene 
and  the  SWNT  surface  is  the  largest.  The  distances  for  the 
benzene  derivatives  are  all  shorter  than  the  distance  of 
3.12  A  for  benzene  (Table  I).  Also,  the  adsorption  energies 
for  the  benzene  derivatives  studied  here  are  larger  than  the 
adsorption  energy  for  benzene,  with  Eads  for  toluene  just 
slightly  above  that  for  benzene  (Table  I),  and  they  are  also 
larger  than  the  Eads  for  the  functional  groups  (when  spin- 
polarization  effects  are  included  for  N02).  All  adsorption 
energies  for  these  aromatics  are  in  the  physisorption  range. 
To  determine  the  relative  contribution  to  the  adsorption  of 
benzene  and  its  functional  groups,  we  now  look  in  detail  at 
the  electronic  structures  of  the  different  systems. 


B.  Electronic  structure 

Here  we  give  results  for  the  electronic  structure  of  the 
benzene  derivatives  as  well  as  of  their  building  blocks.  Fig¬ 
ure  4  shows  the  total  electronic  density  of  states  (DOS)  for 
the  benzene  derivatives  studied  here  and  for  the  small 
closed-shell  molecules  that  contain  the  same  functional 
groups.  The  DOS  for  the  pristine  semiconducting 
SWNT(8,0)  is  given  as  a  reference  in  Fig.  4(b). 

First  we  consider  the  adsorption  of  benzene,  which  is  the 
simplest  aromatic  structure.  Benzene  has  enhanced  stability 
attributed  to  its  delocalized  7 r  electronic  structure  and  its 
closed-shell  structure.  It  has  very  little  effect  on  the  total 
density  of  states  near  the  SWNT  band  gap  [Fig.  4(a)].  The 
molecule’s  highest  occupied  molecular  orbital  (HOMO)  is 
deep  in  the  nanotube  valence  band,  and  its  lowest  uno- 
coupied  molecular  orbital  (LUMO)  is  high  in  the  conduction 
band.  Only  the  states  for  benzene  that  are  about  2  eV  below 
the  gap  are  somewhat  perturbed,  so  the  band  gap  of  the 
SWNT  changed  little.  Mulliken  analysis  shows  a  very  small 
charge  rearrangement  in  the  region  between  the  molecule 
and  the  nanotube  (<0.01  e ),  which  is  within  the  accuracy 
limit  of  the  method  used. 

Although  C6H5N02  is  a  charge  acceptor  due  to  the  elec¬ 
tronegativity  of  N02,  no  charge  transfer  is  found  for  the 


155415-5 


WOODS,  BADESCU,  AND  REINECKE 


PHYSICAL  REVIEW  B  75,  155415  (2007) 


(b)  NT(8,0)+NO2 ,  spin  "down" 


-2.0 


FIG.  5.  (Color  online)  Total 
density  of  states  from  spin- 
polarized  calculations  for  (a) 
NO2/SWNT(8,0),  spin  up;  (b) 
NO2/SWNT(8,0),  spin  down;  (c) 
CH3/SWNT(8 ,0),  spin  up;  (d) 
CH3/SWNT(8,0),  spin  down;  (e) 
NH2/SWNT(8,0),  spin  up;  and 
(f)  NH2/SWNT(8 ,0),  spin  up. 


nitrobenzene/SWNT.  The  HOMO  of  nitrobenzene  is  deeper 
in  energy  than  the  HOMO  of  N02,  and  it  is  completely  filled. 
Figure  4(c)  for  C6H5N02/SWNT  shows  a  perturbation  of  the 
total  DOS  at  about  1.5  eV  above  the  band  gap  due  to  the 
LUMO  of  the  molecule.  Otherwise,  there  is  little  effect 
around  the  band  gap  and  Fermi  level  from  the  N02  func¬ 
tional  group.  If  a  stronger  hybridization  were  to  take  place, 
one  would  expect  that  the  equilibrium  configuration  of  the 
molecule  would  be  closer  to  the  transverse  direction  due  to  a 
stronger  attraction  of  the  N02  to  the  SWNT.  This  was  not 
observed  here. 

For  the  toluene  adsorption,  it  is  evident  from  Fig.  4(e)  that 
the  semiconducting  nature  of  the  SWNT  is  also  preserved 
and  that  the  energy  gap  changed  very  little.  We  notice  that 
the  profile  of  the  total  DOS  is  very  similar  to  that  for 
C6H6/SWNT.  This  is  consistent  with  the  adsorption  energy 
(Table  I)  being  close  to  that  for  benzene  and  to  the  fact  that 
CH3  has  little  electronegativity.  No  charge  transfer  is  found 
between  the  toluene  and  the  SWNT,  and  the  adsorption  re¬ 
sults  primarily  from  the  tt-tt  interaction. 

Aniline  is  expected  to  be  a  charge  donor  because  NH2  is 
electropositive.  Here  a  small  charge  redistribution  is  found 
between  this  molecule  and  the  SWNT.  Figure  4(g)  indicates 
that  there  are  no  other  hybrid  states  in  the  semiconducting 
gap.  However,  a  different  hybrid  state  appears  below  the  top 
of  the  valence  band  but  does  not  affect  the  Fermi  level.  A 
decomposition  of  the  DOS  indicates  that  the  state  is  formed 
due  to  hybridization  between  the  highest  SWNT  valence 
band  and  the  aniline  HOMO  and  that  it  is  mainly  of  p  char¬ 
acter. 

In  order  to  understand  better  the  adsorption  of  the  ben¬ 
zene  derivatives  on  nanotubes,  we  now  analyze  the  elec¬ 
tronic  structure  of  the  different  functional  groups  themselves. 
We  recognize  that  the  unsaturated  functional  groups  have  a 
qualitatively  different  electronic  shell  structure  than  when 
they  are  saturated  by  attaching  to  C6H5  (from  the  benzene 


ring)  or  other  fragments.  The  functional  groups  N02,  NH2, 
and  CH3  considered  here  all  have  a  partially  filled  valence 
electron  shell.  When  they  are  substituted  for  a  H  atom  in 
benzene,  the  resulting  complex  has  a  filled  valence  shell. 
Thus  the  functional  groups  will  adsorb  differently  on  the 
nanotube  surface  if  they  are  attached  to  a  benzene  ring  than 
when  they  are  unsaturated.  To  clarify  this  issue,  we  studied 
the  adsorption  of  simple  closed-shell  molecules  that  contain 
the  functional  groups.  Then  we  describe  the  open-shell  frag¬ 
ments  themselves. 

In  Figs.  4(d),  4(f),  and  4(h),  the  total  DOS  is  given  for 
nitromethane  CH3N02,  methane  CH4,  and  amonia  NH3.  It 
can  be  seen  that  the  adsorption  of  these  molecules  also  has 
little  effect  on  the  density  of  states  of  the  clean  nanotube 
around  the  SWNT  band  gap.  The  compounds  have  a  HOMO 
and  a  LUMO  far  in  the  valence  and  conduction  bands  of  the 
SWNT(8,0),  respectively.  In  addition,  their  closed-shell 
structure  makes  them  unreactive  to  other  compounds  in  a 
similar  manner  as  the  closed-shell  benzene  derivatives.  In 
summary,  we  found  that  all  closed  valence  shell  adsorbates 
examined  here  adsorb  weakly  and  do  not  change  the  nano¬ 
tube  electronic  structure  significantly. 

The  total  density  of  states  for  the  partially  filled  valence 
shell  fragments  NH2,  N02,  and  CH3  are  given  in  Fig.  5. 
Their  behavior  is  different  from  the  closed-shell  systems.  For 
these  systems,  spin-polarized  calculations  have  been  done. 
The  HOMO  of  N02  is  close  to  the  valence-band  edge  of  the 
SWNT  and  is  partially  filled.  After  the  adsorption,  this  state 
is  hybridized  with  the  highest  conduction  states  of  the  nano¬ 
tube  and  stays  under  the  band  gap.  Also  a  peak  appears  right 
above  the  Fermi  level  for  “spin  up.”  This  state  is  empty  and 
contains  little  admixture  from  the  C  orbitals.  No  such  state  is 
found  for  “spin  down.”  The  Mulliken  analysis  shows  that 
there  is  very  little  charge  transfer  in  the  system  —0.01  e.  For 
NH2  we  also  found  a  peak  mainly  of  p  character  from  the 
fragment  400  meV  above  the  Fermi  level  for  spin  up,  which 
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FIG.  6.  Plots  of  charge  density  in  five  subbands  at  the  top  of  the 
valence  bands  for  [(a)  and  (b)]  C6H5NH2/SWNT(8,0)  and  [(c)  and 
(d)]  C6H5NO2/SWNT(8,0).  In  each  figure  the  dark  areas  represent 
the  highest  densities.  The  contour  plots  are  on  a  logarithmic  scale, 
(a)  and  (c)  show  plots  for  the  SWNT  axial  cuts;  (b)  and  (d)  show 
plots  for  cuts  perpendicular  to  the  SWNT  axis  through  the  func¬ 
tional  atoms  H  for  aniline  and  O  for  nitrobenzene. 

is  also  an  empty  state.  No  such  peak  is  obtained  for  spin 
down. 

For  CH3,  the  spin-polarization  calculations  result  in 
spin-up  states  right  under  the  Fermi  level  that  are  completely 
filled  and  in  spin-down  states  that  are  empty.  These  states 
result  from  a  strong  hybridization  between  the  CH3  and  the 
nanotube.  This  bond  is  of  sp 3  character  and  is  a  result  of  the 
admixture  between  the  C  states  from  the  radical  and  the  C 
states  from  the  four  atoms  from  the  nanotube  under  the  CH3. 
We  also  found  that  there  is  a  deformation  of  the  nanotube 
surface  under  the  radical  upon  its  adsorption,  which  is  a 
manifestation  of  the  sp 3  bonding.  Mulliken  analysis  esti¬ 
mates  that  there  is  a  charge  transfer  of  0.2  e  toward  the  mol¬ 
ecule.  For  the  case  of  CH3,  we  obtain  stronger  bonding,  more 
hybridization,  and  more  nanotube  deformation  than  for  the 
other  adsorbates  studied  here.  This  results  from  the  strong 
C-C  bond  formed  with  the  carbon  nanotube. 

Further  insight  can  be  gained  from  charge-density  contour 
plots  of  the  electronic  densities.  They  were  obtained  with  the 
LEV00  computer  code.44  In  Fig.  6  we  present  contour  plots 
for  C6H5N02  and  C6H5NH2.  For  nitrobenzene/SWNT,  the 
charge  density  is  simply  a  superposition  of  the  charge  densi¬ 
ties  for  the  molecule  and  the  SWNT  separately,  and  very 
little  hybridization  occurs.  Nitrobenzene  sits  on  the  top  site 
and  interferes  little  with  the  SWNT  electronic  structure.  For 
aniline/SWNT,  hybridization  occurs  between  the  highest  va¬ 
lence  bands  just  below  the  SWNT  gap  and  the  aniline 
HOMO.  The  n t  bonds  of  the  C  atoms  from  the  nanotube  are 
somewhat  distorted  and  are  hybridized  with  the  N  levels 
from  aniline.  The  H  atoms  from  the  NH2  functional  group 
from  aniline  do  not  affect  the  SWNT;  the  hybridization  oc¬ 
curs  mainly  with  the  N  atom. 

Figure  7  shows  total  charge-density  plots  for  the  axial  and 
transverse  nanotube  directions  for  the  spin-polarized  calcula- 


FIG.  7.  Total  charge-density  plots  for  [(a)  and  (b)] 
NO2/SWNT(8,0)  and  [(c)  and  (d)]  CH3/SWNT(8 ,0).  The  contour 
plots  are  on  a  logarithmic  scale.  In  each  figure  the  areas  with  the 
largest  density  are  represented  by  dark  gray,  (a)  and  (c)  show  plots 
for  the  SWNT  axial  cuts;  (b)  and  (d)  show  plots  for  cuts  perpen¬ 
dicular  to  the  SWNT  axis  through  the  functional  atoms  N  and  C. 


tions  for  N02/SWNT  and  CH3/SWNT.  For  the  N02  case,  the 
states  above  the  Fermi  level  have  very  little  admixture  from 
the  nanotube,  and  thus  no  charge  transfer  between  the  N02 
and  nanotube  orbitals  occurs.  For  CH3,  strong  sp3  hybridiza¬ 
tion  of  the  C  atom  with  nanotube  states  occurs,  resulting  in  a 
significant  charge  transfer  in  the  system. 

IV.  CONCLUSIONS 

The  present  results  from  the  DFT  calculations  can  be  used 
to  determine  the  important  physical  features  in  the  adsorp¬ 
tion  process  of  benzene  derivatives  on  SWNTs.  One  factor  is 
that  the  noncovalent  tt-tt  interaction  between  the  benzene 
ring  and  the  nanotube  is  the  dominant  adsorption  mechanism 
regardless  of  the  charge-transfer  properties  of  the  functional 
groups  or  the  locations  of  their  electronic  states  with  respect 
to  the  nanotube  gap.  The  functional  groups  have  only  a  sec¬ 
ondary  effect  on  the  adsorption  and  have  small  effects  on  the 
electronic  properties  of  the  semiconducting  SWNT.  The  cal¬ 
culated  binding  energies  and  equilibrium  distances  are  con¬ 
sistent  with  physisorption. 

In  the  case  of  the  adsorption  of  the  open-shell  functional 
fragments  studied  here,  the  existence  of  a  partially  filled  or¬ 
bital  and  the  position  of  the  HOMO-LUMO  with  respect  to 
the  nanotube  gap  are  also  other  important  factors  in  the 
nanotube  adsorption.  We  see  that  for  simple  closed-shell 
molecules  similar  to  the  open-shell  fragments,  the  adsorption 
is  similar  to  that  for  the  benzene  derivatives.  On  the  other 
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hand,  the  adsorption  of  the  open-shell  fragments  themselves 
can  change  the  nanotube  electronic  structure  significantly. 

The  comparison  of  the  equilibrium  distances  and  energies 
for  the  different  adsorbed  aromatic  molecules  suggests  some 
competition  between  the  benzene  ring  and  the  functional 
groups.  The  equilibrium  distances  for  the  benzene  deriva¬ 
tives  are  smaller  than  those  for  C6H6  and  greater  than  those 
for  the  respective  functional  groups,  with  the  exception  of 
toluene  (Table  I).  Also,  the  adsorption  energies  for  the 
benzene-derived  aromatics  are  larger  than  the  adsorption  en¬ 
ergy  for  benzene,  indicating  that  functional  groups  affect  the 
adsorption  process  to  some  extent.  Again,  the  exception  is 
toluene,  for  which  the  energy  is  very  close  to  that  for  ben¬ 
zene  (Table  I).  The  carbon  ring  of  benzene  tends  to  keep  the 
molecule  at  a  distance  characteristic  of  tt-tt  stacking  in  gra¬ 
phitic  structures,  whereas  the  functional  group  tends  to  bring 
the  molecule  closer  to  the  SWNT. 

These  results  suggest  that  these  benzene  derivatives  can 
be  used  for  the  noncovalent  functionalization  of  carbon 


nanotubes.  These  aromatics  might  also  play  a  role  in  systems 
for  biological  species  immobilization  on  nanotubes.  Our  re¬ 
sults  for  the  adsorption  of  open-shell  systems  might  be  im¬ 
portant  for  the  chemical  functionalization  of  carbon  nano¬ 
tubes. 
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